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ABSTRACT 

I report on recent progress in the study of mass generation in non-abelian gauge 
theories, and concentrate on the origin of constituent masses in QCD. I argue that 
a consistent formalism has been developed to describe chiral symmetry breaking via 
multiplicatively renormalisable Dyson-Schwinger equations. Its main consequence 
is the existence of a critical value of the coupling beyond which massless solutions 
do not exist. 



It is well-known that quarks are described by two different kinds of masses. First 
of all, there are the masses which enter as parameters in the lagrangian, or in sum- 
rules, the current masses. For light quarks, these are small parameters (m„ ~ 0.005 
GeV, m<i ~ 0.01 GeV). On the other hand, the kinematics of quarks {e.g. in B 
or D decays) and the masses of hadrons imply much larger masses for quarks, the 
constituent masses, of the order of 0.3 GeV for u and d quarks. These masses come 
from the interactions of quarks propagating in a non-trivial vacuum, similar to the 
interactions of electrons with the background held of a crystal. Hence 97% of the 
observed hadronic mass comes from the interaction of quarks with the QCD vacuum. 

In the limit where one neglects current masses, the SU{3) symmetry of the QCD 
lagrangian gets doubled into a left- and a right- SU (3), as the left-handed and right- 
handed chiralities do not mix without a mass term. Chirality then becomes a con- 
served quantity, and the QCD lagrangian acquires a chiral symmetry. One then 
expects the dynamics of the theory to break this symmetry into and SU{3) with 
mass terms, appearing as poles in the quark propagators. We shall write the quark 
propagator as 
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Note that we have Wick-rotated to Euclidean space, where we shall work hereafter. 
Chiral symmetry breaking occurs if E(& 2 ), or equivalently G(k 2 ) } are nonzero. Note 
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that the question we are addressing is that of mass generation in gauge theories, and 
the formalism we are going to develop can be applied to any gauge theory (e.g. to 
technicolour). 

Let us hrst explain why Dyson-Schwinger equations are relevant in this context. 
First of all, perturbation theory, at least in lowest order, fails to produce a mass 
for massless objects: the one-loop corrections to the propagator merely multiply 
the propagator I/7 • p by a function of p 2 } and do not shift the position of the 
pole. In fact, a renormalisation group argument probably due to Callan, Dashen 
and Gross [1], indicates that dynamical mass generation lies in the non-perturbative 



regime of QCD. Indeed, the quark mass m = wE(m 2 ) can only be a function of the 
coupling constant and of the renormalisation point //: m(g } fi) = fiM(g) } with M a 
dimensionless function. If the mass is to be physical, it cannot depend on //, hence 
fidm/dfi = m + (3(g)dm/dg = 0, with (3 = dg/dfi. In QCD with two flavours, the 
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one-loop (3- function is (3(g) ~ — t^t*/ 3 , which leads to: 
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Such a dependence on as cannot be produced via perturbation theory, even if re- 
summed, as m does not have an asymptotic expansion as as — > 0. Hence the problem 
is intrinsically non-perturbative. A caveat is here in order. It was mentioned above 
that perturbation theory fails to produce a mass. This is true only if one is careful 
about gauge invariance. Any additional term, regularization method, or resummation 
ansatz that breaks gauge invariance will in general produce a mass term, which will 
come, not from the physics, but rather from the clumsiness of the method. Hence the 
preservation of gauge invariance is of utmost importance in this problem. 

Field theory contains few nonperturbative statements. Some of these, particularly 
suited for the present study, are the equations of motion of Green's functions, the 
Dyson-Schwinger (DS) equations. We sketch in Fig. 1 the equation for the two-point 
function, S(p 2 ) which we expect to have a pole at the physical mass of the quark. 




Figure 1: A pictorial representation of the Dyson-Schwinger equation for quarks. 
The exact two-point functions are indicated by a hatched circle, the exact three- 
point function by a doubly hatched circle. 

Note that we work in the axial gauge, hence no ghost diagrams are present. This DS 
equation relates the exact quark propagator to two a priori unknown quantities: the 
exact gluon propagator and the exact three-point function. We shall assume that the 
gluon propagator is independently known, e.g. from the gluon DS equation in the 



quenched approximation [2], and we shall write it as 

DM) = <l 2 D(q 2 )Dl( q 2 ) (3) 

with D° the perturbative gluon propagator. We shall here rather concentrate on the 
approximations one can make to the three-point function. One has to guess its form, 
and three levels of approximation are possible. First of all, the simplest assumption 
is to take the perturbative vertex and assume that this is good enough. This leads 
to the "rainbow approximation", and to the breakdown of gauge invariance. The 
second level of approximation [3] is to use another non-perturbative statement, the 
Ward-Takahashi (WT) identities, and choose a vertex that will respect these, and 
hence not explicitly break gauge invariance. As we shall see, this is not restrictive 
enough, and one needs to go to the third level in the imposition of gauge invariance, 
and demand one further constraint [4], namely the matching to perturbation theory 
at high q 2 } in order to obtain a sensible equation for chiral symmetry breaking. 

Let us hrst concentrate on WT identities, as shown in Fig. 2. One can see that 
the projection of the three-point function onto the gluon momentum is related to the 
difference of two inverse propagators. Hence the part of the vertex function which is 
parallel to the gluon momentum, i.e. the longitudinal part, is known. The choice of 
vertex function is then reduced to guessing the transverse part IV, which does not 
enter the WT identities. 
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Figure 2: A pictorial representation of the Ward-Takahashi identities using the same 
convention as in Fig. 1. 

The simplest assumption [3] is to take Tt = 0. One then obtains [5] two linear 
decoupled equations for F and G: 

1 = F{q 2 ) + a s F{q 2 )jdk 2 K 1 {k\q 2 ) + asjdk 2 K 2 {k\q 2 )F{k 2 ) (4) 

= G(q 2 ) + a s G(q 2 ) [ dk 2 K^k 2 ,q 2 ) + a s I dk 2 K 2 (k 2 ,q 2 )G(k 2 ) (5) 



The kernels Ki and K 2 are the same in each equation, and can be found in Ref. [5]. 
They are both proportional to the exact gluon propagator, and hence to D(q 2 ) of 
Eq. (3). Their most important property is that the integral of Ki diverges loga- 
rithmically for large k 2 : one encounters divergences similar to those of perturbation 
theory, which need to be renormalised away. 



The renormalisation proceeds along the lines of the usual prescription: one intro- 
duces a renormalisation point // such that F(q 2 ) = Zj(fi 2 )Fji(q 2 ) } and one chooses 
fj 2 Fr(jj 2 ) = 1. As there is no quark mass term in our lagrangian, the mass squared E 
of Eq. (1) does not get renormalised. This amounts to saying that G and F have the 
same renormalisation factor Zj(/j, 2 ). One also has to renormalise the gluons, which 
is done by similarly renormalising the gluon propagator (3): D(q 2 ) = Z(jj 2 )Dr(jj 2 ). 
The (infinite) renormalisation constant Z(jj 2 ) is absorbed in the definition of the 
renormalised coupling. One then subtracts from Eqs. (4,5) their value at q 2 = fj 2 and 
this subtraction removes the remaining renormalisation constant Zf(jj 2 ). 

One runs however into several problems: 

• In the usual multiplicative renormalisation (MR), one defines the renormalised 
coupling as c^(ji) = Z(fj,)a b g re , and this leads to the renormalisation group. Here 
however, this definition is not sufficient to remove Z(fi) from the DS equation. One 
has to use a definition that breaks MR: 

R _ Z{ji)a h i re 
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As the kernel K\ is specific to the quark equation, one looses the universality of the 
renormalised QCD coupling: gluons couple differently from quarks! Furthermore, one 
also obtains the wrong asymptotic behaviour for as(q 2 ) at large q 2 . 

• The renormalised equation for G(q 2 ) at large q 2 is dominated by the terms resulting 
from Ki, and one can show that these behave like log(log(g 2 )). Hence the leading 
terms of Eq. (5) look like G(q 2 ) = as(fi)G(q 2 )log(log(q 2 )) } and the equation has no 
consistent nonzero solution. 

• If one proceeds and solves the equation for F, one obtains a solution for any as(fi) 
and these solutions have a pole at the origin. Hence the quarks are massless and 
unconhned. 

All these problems come from the large k 2 region, which is where we expect 
perturbation theory to hold. In the region where q 2 ss k 2 >>> p 2 } one can calculate 
the one-loop perturbative corrections to the vertex, and hence, by subtracting the 
longitudinal vertex derived from the WT identities, get an asymptotic expression for 
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where £ is the quark anomalous dimension. Furthermore, one has the constraint that 
by definition, for any p and q } Tt • (q — p) = 0. One can then show that the general 
form of IV, for E = 0, is: 

S-^ T S = [q 2 F(q 2 ) - /)F(/)] 7At(g2 -^ ) + (gAt +^ )(7 ^- 7 - g) (8) 

where T> behaves like q 4 at large q 2 } is symmetric in p and q and must be nonsingular. 
We choose in the following T> = (q 2 + p 2 ) 2 . 



We now concentrate on the chiral solution G = 0. The equation is similar to 
Eq. 4, with new kernels K\ and A 2 replacing K\ and A 2 . These kernels are again 
given in Ref. [5]. Their essential property is that the divergence that used to be in the 
K\ term is now shifted to A 2 . As a consequence, the definition of the renormalised 
as, Eq. (6), now contains a finite integral over Ai, and hence differs from the usual 
MR prescription by a finite calculable factor, hence the usual redefinition a^(ji) = 
Z(ji)a h g re now works and produces an equation free of singularities. 

The renormalised equation looks as follows: 



q 2 )q 2 F R (q 2 ) = <f>(» 2 ) + dk 2 [K 2 (k 2 ,q 2 ) - K 2 (k 2 , » 2 )]F R (k 2 ) 



(9) 



with 
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This equation has two regimes. If (f>(q 2 ) > for all q 2 } then there is a solution, which 
means that F ^ 0, G = is allowed and that chiral symmetry remains unbroken. On 
the other hand, if (f>(q 2 ) = for some value of q 2 then there is no solution for F and 
G = is not allowed. 

To study the positivity of <f> for a general gluon propagator, one can use a Kallen- 



Lehmann representation for the gluon propagator, D(k 2 
fact that the kernel K\ is linear in D(k 2 ) } Ki(q 2 } k 2 
rewrite 
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H(cr } q 2 ) > implies that one has a massless solution. We show in Fig. 3 the region in 
(as, a) space for which the function H is positive definite, and the upper curve shows 
the values of as(jJ 2 ) for which H(cr } fi 2 ) = 0, i.e. for which the effective coupling (6) 
becomes infinite, and where one expects the equation not to have a solution. 
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Figure 3: The shaded region is where the function H is positive definite for all q ; 
the thick curve shows the values of as for which H[<7, fi ) = 0. 

Hence we see that massless solutions exist only for small values of cts(fi), and that 
there is a critical value of as(fi) (one should really speak of a renormalisation-group 



invariant such as cts(fi)fi 2 Dr(/j, 2 )) beyond which one ceases to have massless solutions. 
Futhermore, the chiral symmetry of QCD can also be destroyed by the presence of 
massive, high-cr modes in the gluon propagator. Although the exact value of a c g tUcal 
will depend on the gluon Kallen-Lehmann density p(cr) } it is clear that a critical value 
of the coupling will always exist. 

To show what happens to the massless solution in an explicit case, we can as a 
toy-model use the gluon propagator derived in Ref. [2]: 
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The equation leading to this solution followed the old BBZ approach [3] and hence 
violated MR, and did not give rise to the renormalisation group. As a result, we were 
able to derive a solution only for a given cts(fi), «s(/i) = 1.4, as the renormalisation 
group which compensates a variation of as by a variation of D(k 2 ) did not work. The 
solution has a cut behaviour: near k 2 = 0, D(k 2 ) ~ (A; 2 ) -0 ' 22 , and hence corresponds 
to confined gluons, with no pole near the origin. As one expects the pole to be tapered 
by nonperturbative effects, and as a theorem [3] requires the axial gauge propagator 
to be infinite at the origin, one expects an improved solution satisfying MR to look 
roughly the same. Hence we use the solution found in [2] and the best we can do 
at present is to change the value of as(fi) only in the quark equation. One then 
obtains the curves of Fig. 4. They show that as the value of as(fi) increases, the 
quark propagator experiences oscillations, until eventually no solution can be found, 
for a critical value of as of the order of 1.5. It is interesting that the value of as that 
we obtained from the gluon equation is very close to the critical value. 
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Figure 4: The solution for massless quark propagators that come from the non- 
perturbative gluon propagator of Ref.l, for as(fi) =0.2 (plain), 0.6 (dashed), 1.0 
(dot-dashed) and 1.4 (dashed). 

To conclude, we now have a formalism to handle the Dyson-Schwinger equations 
which takes into account everything we know about QCD: 

• The renormalisation is multiplicative, and hence gives rise to the renormalisation 
group equations; 

• Gauge invariance is enforced via the WT identities, and via the imposition of a 



transverse vertex that matches perturbative QCD in the ultraviolet region; 
•The quark propagator matches smoothly with the perturbative ansatz at large mo- 
mentum. 

This formalism predicts that one keeps massless solutions only for small as, and 
hence there is a critical value of as beyond which a mass must be generated for the 
equation to have a solution. One can also show that the equation for the mass term 
has a consistent behaviour at large momentum. 

One of course now needs to find the solutions in the massive phase and to extend 
the present formalism to the gluon case. One will then have a consistent held theoretic 
description of chiral symmetry breaking at the quark level. 
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